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Abstract: 

A generalization of the AdS/CFT conjecture postulates a duality between IIA string theory 
and 16 supercharge Yang-Mills quantum mechanics in the large N 't Hooft limit. At low 
temperatures string theory describes black holes, whose thermodynamics may hence be stud- 
ied using the dual quantum mechanics. This quantum mechanics is strongly coupled which 
motivates the use of lattice techniques. We argue that, contrary to expectation, the theory 
when discretized naively will nevertheless recover continuum supersymmetry as the lattice 
spacing is sent to zero. We test these ideas by studying the 4 supercharge version of this 
Yang- Mills quantum mechanics in the 't Hooft limit. We use both a naive lattice action and 
a manifestly supersymmetric action. Using Monte Carlo methods we simulate the Euclidean 
theories, and study the lattice continuum limit, for both thermal and non-thermal periodic 
boundary conditions, confirming continuum supersymmetry is recovered for the naive action 
when appropriate. We obtain results for the thermal system with N up to 12. These favor 
the existence of a single deconfined phase for all non-zero temperatures. These results are an 
encouraging indication that the 16 supercharge theory is within reach using similar methods 
and resources. 
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1. Introduction 



Despite progress in understanding quantum gravity there has been relatively little progress in 
theoretically realizing a quantum black hole, the key object in quantum gravity. Substantial 
progress was afforded by the work of Strominger and Vafa [1] who showed that for certain 
black holes in string theory the constituent microstates can be explicitly counted and the 
result is in agreement with the usual Bekenstein-Hawking entropy [2]. The calculation relies 
on a continuation to a weakly coupled theory where states can be explicitly counted, and 
this number is then argued to be independent of the continuation. This gives an indirect 
understanding of the microstates of the black hole. Another beautiful result of Strominger 
is that the entropy of black holes with an AdSs near horizon geometry can be computed 
on very general grounds as the gravitational theory describing the black hole can be argued 
to be a 2-d conformal field theory [3]. One does not know the theory, but does know the 
central charge, which is sufficient to compute the entropy using Cardy's result. Thus again 
we gain insight into the microstates, although we are not able to solve the theory governing 
the quantum black hole. For a review of these approaches we refer the reader to [4] and for 
an interesting new proposal concerning understanding the microscopic description without 
resorting to a weak coupling continuation see the work of Mathur [5] . 

A very promising new direction for the study of quantum gravity is the AdS/CFT corre- 
spondence and its generalizations [6, 7]. Here one can describe certain black holes in terms of 
the worldvolume theory of the D-branes that compose them. These are 16 supercharge gauge 
theories in various dimensions, taken in the large N 't Hooft limit and at finite temperature. 
The regime in which they describe string theory black holes is one in which they are strongly 
coupled, and solving this theory would allow one to directly study the quantum properties 
of the dual black hole, including its thermodynamic properties. The problem of computing 
directly in this strongly coupled theory is a technical one, not one of principle. While analytic 
methods have made much progress in describing the planar limit of the 4 dimensional version 
of this correspondence using integrability techniques [8] it is not generally understood how 
to go beyond this limit. It particular black holes are precisely objects that live beyond this 
limit, having dual energy densities of order 0(N 2 ). Even protected quantities that can give 
information about the strongly coupled theories through weak coupling calculations have not 
yet revealed information about black holes [9] and if they do in the future, they will remain 
indirect probes of quantum black hole physics. Another direction for direct calculation is 
the use of the Gaussian approximation, pioneered in this context for Yang-Mills quantum 
mechanics in [10]. Certain agreement with black hole thermodynamics was extracted [11], 
although the approximation scheme is essentially uncontrolled, having no small parameter. 

The aim of this paper is to investigate the feasibility of using lattice methods to directly 
simulate the strongly coupled worldvolume theories. Hiller et al [12] have used lightcone meth- 
ods to give evidence that the graviton propagator is correctly reproduced in the 2 dimensional 
version of the correspondence. Here, we test the use of direct Monte Carlo simulation of the 
lattice regulated path integral to investigate these systems. Since we are interested in study- 
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ing the theory at finite temperature, Monte Carlo techniques seem the most natural as they 
allow us to study observables directly in the thermal ensemble. In addition, Monte Carlo 
methods are often the only feasible approach in higher dimensions. 

These theories are necessarily supersymmetric, and this presents technical challenges for 
the lattice. The first problem encountered in such a study is how best to discretize such a 
supersymmetric theory. There has been much interest in the lattice study of supersymmetric 
theories. Aside from our current application, supersymmetry plays a prominent role in efforts 
to construct theories which go beyond the standard model of particle physics and the lattice 
provides the possibility of studying nonperturbative effects such as dynamical supersymmetry 
breaking, of relevance for phenomenology. Unfortunately, in many cases, supersymmetry is 
broken at the classical level in such discretizations and can only be regained in the quantum 
continuum limit after a great deal of fine tuning in the lattice theory. Recently significant effort 
has gone into attempts to construct lattice theories which retain an exact supersymmetry 
at non-zero lattice spacing. Two approaches have been used; lattice models arising from 
orbifolding a supersymmetric matrix model [13, 14, 15, 16, 17] and constructions based on 
discretizing a topological or twisted form of the continuum supersymmetric theory [18, 19, 
20, 21, 22]. In the case of Yang-Mills theories it appears that these approaches are intimately 
connected [23, 24]. The approach based on twisting has also been studied in the case of 
Wess-Zumino and sigma models [25, 26, 27, 28, 29]. The philosophy behind these approaches 
is that this residual lattice supersymmetry will help to protect the theory from the dangerous 
radiative corrections which lead to fine tuning. 

However, in a sufficiently low number of dimensions the super renormalizable nature of the 
models leads to much reduced fine tuning problems - typically there are only a finite number 
of divergences in such theories which may be computed in perturbation theory [30] . In the case 
considered here, supersymmetric Yang-Mills quantum mechanics, we argue that no relevant 
supersymmetry breaking counterterms can be written down. Naive discretizations of such a 
theory, while breaking supersymmetry classically by terms of the order of the lattice spacing, 
flow automatically to the supersymmetric theory in the naive continuum limit corresponding 
to a vanishing lattice spacing. 

In this paper we have examined the effectiveness of using Monte Carlo simulation to 
study super Yang-Mills quantum mechanics in the large iV 't Hooft limit in the case of the 
4 supercharge model. This theory, while exhibiting many features which are qualitatively 
similar to its 16 supercharge cousin, is computationally easier and a useful 'warm-up' exercise 
before tackling the 16 supercharge case. It has been studied using Hamiltonian methods 
in [31, 32] for SU{2). We have studied two different discretizations of this theory. One of 
these possesses 2 exact lattice supersymmetries corresponding to the dimensional reduction 
of the twisted model derived in [21]. The other corresponds to a naive discretization of the 
continuum theory constrained only by the necessity of employing a lattice derivative which 
forbids fermion doubling. Both are invariant under lattice gauge transformations. 

We give analytical arguments that the naive action should not require fine tuning to 
attain the correct continuum limit. We reinforce this conclusion by conducting numerical 
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simulations of the model on finite lattices at zero temperature. The results of these simulations 
show that the vacuum energy approaches zero as the lattice spacing is reduced and that 
expectation values computed in this naive lattice theory agree with those computed with 
the supersymmetric lattice action. Since the simulations of the supersymmetric system are 
much more computationally demanding we have primarily used the naive action in subsequent 
simulations of the thermal system. 

To set the stage for the interpretation of our thermal results we have also simulated 
the quenched theory. These confirm the previous expectation of a rapid crossover in the 
thermodynamic behavior of the system at some temperature which becomes a sharp phase 
transition as N — > oo [33, 34]. For high temperatures the theory appears deconfined but 
below the critical temperature the system enters a confining phase with non-zero vacuum 
energy. Our new results from simulations of the full 4 supercharge theory however favor a 
simple phase structure with a single deconfined phase for all non-zero temperatures - as we 
would expect for the behaviour of the 16 supercharge theory from holographic considerations. 
We are able to simulate this 4 supercharge thermal theory with iV up to 12, already sufficient 
to see the asymptotic 't Hooft scaling. 

The plan of the paper is as follows. In section 2 we introduce the Yang-Mills model in 
the continuum including a discussion of its renormaliz ability and vacuum structure. Section 3 
describes the two lattice actions we have utilized while we present our results for the quenched, 
zero temperature and thermal theories in section 4. Our final section discusses the implication 
of our results both in the context of the 4 supercharge model and its string theory cousin 
with 16 supercharges. 

Finally, the paper ends with several appendices containing more technical details on 
holographic duality, the quantum corrections to the classical moduli space of the theory, 
and further numerical results concerning the extrapolation of our finite lattice data to the 
continuum limit. It also contains a detailed description of the RHMC algorithm used in our 
simulations. 

While this paper was in preparation we received a paper which utilizes momentum space 
methods to study the same system for S'C/(4) with results that are in at least qualitative 
agreement with ours [35]. 

2. 4 supercharge Yang-Mills quantum mechanics 

The Euclidean supersymmetric SU(N) Yang-Mills quantum mechanics we are interested in 
can be thought of as arising from the classical dimensional reduction of N = 1 super Yang- 
Mills in 4 dimensions. The matter content in the quantum mechanics arises from reduction 
of the gauge field in 4-d, giving 3 adjoint bosonic scalar fields Xi, where i = 1, ... ,3 and 
the adjoint gauge field, A. The 4-d adjoint fermions can be equivalently written in either 
a Weyl or Majorana representation. Here we will use the Weyl representation, giving a 
complex 2-component fermion field, Vl/ Q with a Euclidean spinor index a, transforming in the 
adjoint of the gauge group. Dimensionally reducing this, the spinor index becomes an internal 
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symmetry. The quantum mechanics Euclidean action is given as, 

5 = ^NTij R dT \^{D T X t ) 2 - l - [X^Xj] 2 + m r D T V - W [Xi,*]\ , (2.1) 

with r being compact Euclidean time with radius R. The covariant derivative is defined as 
D T = d T + % L4(t), •] on an adjoint field, and we use the usual (Euclidean) chiral conventions, 

with a 1 being the Pauli matrices. We have written the SU(N) gauge theory using the 't 
Hooft coupling A = Ngy M , with gyu the Yang-Mills gauge coupling. The supersymmetry 
transformation derives from that of the reduction of the parent 4-d = 1 super Yang-Mills, 
giving, 

SA = -i^a T £ 
5Xi = -Ma^ 

8* = 2 (<t™£) (D T Xi) + i (a^'e) [X h X,] (2.3) 

where we use the notation a^ v = \ (a^a u — a^a^ 1 ). The theory has a global SO (3) symmetry 
which can be interpreted as a rotational invariance in the target space of the theory, with 
action, 

X 'i = A (i)A K = Af 1)Q *^ (2.4) 

where A^) gives the representation of 50(3) acting on vectors, and A(i/2) the representation 
on spinors. 

Since the Euclidean time is topologically 5 1 we have two spin structures which give either 
periodic or anti-periodic boundary conditions for the fermion ^ as we traverse the time circle. 
In the latter case, the path integral Z is the usual thermal partition function with temperature 
T, 

Z a {R) = Tre- RH (2.5) 

where R is the inverse temperature 1/T, and H is the Hamiltonian operator. Taking periodic 
boundary conditions, 

Z p {R) = Tr{-l) p e- Rfl (2.6) 

where F is the fermion operator, and hence Z p gives the Witten index, whose value counts 
the difference in the number of fermionic and bosonic ground states H\ip >= 0, and hence 
should not depend continuously on R. 

2.1 The 16 supercharge theory and its IIA closed string theory dual 

The 16 supercharge theory, which arises as the 1 dimensional reduction of M = 1 super Yang- 
Mills in 10 dimensions, takes a very similar form to the above. The only difference is that 
the reduction of the 10-d gauge field gives rise to 9 adjoint scalars Xj, with i = 1, . . . , 9, and 
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the 10-d Majorana fermions yield an adjoint fermion ^ now with spinor index p, = 1, . . . , 16. 
Note that unlike the 4 supercharge theory, there is no Weyl representation available for the 
fermions, and hence one obtains a Pfaffian from integrating out these fields rather than a 
determinant. 

Since the 16 supercharge theory is so similar to the 4 supercharge theory, one might ask 
why we have decided to simulate the 4 supercharge case here. Clearly the 16 supercharge 
case will generate a fermion operator which is four times larger than its 4 supercharge cousin. 
With current simulation algorithms for near massless dynamical fermions this will lead to 
a factor of 16 slowdown in our simulations. In addition even in the continuum Euclidean 
theory the Pfaffian that results after integration over the fermions is in general complex as 
seen from considering the zero momentum sector of the theory [36] . In contrast to this the 
continuum 4 supercharge theory has a Pfaffian which can be recast as a positive definite real 
determinant, and we later show that our naive lattice discretization also has this property. 
The feasibility of using Monte Carlo methods to simulate such a system with complex Pfaffian 
depends upon whether the resulting phase fluctuations are relatively small and infrequent as 
we approach the continuum limit. In the latter case standard reweighting techniques can be 
used to evaluate expectation values generated in the phase quenched ensemble (see eg. [37]. 
Because of these technical subtleties we have decided to initially study the 4 supercharge 
theory before moving to the 16 supercharge case. 

Following the AdS-CFT conjecture Itzhaki et al [7] have argued that SU(N) 16 su- 
percharge Yang-Mills theories taken in the large N 't Hooft limit are dual to certain closed 
superstring theories in the near horizon region of N coincident D-branes. In particular SU(N) 
super Yang-Mills quantum mechanics is supposedly dual to the IIA string theory describing 
N DO-branes. We review this correspondence in more detail in the appendix A. Here we 
simply state the results that the analysis gives. 

At a finite temperature T = 1/R the quantum mechanics theory then describes a gas of 
N DO-branes in the dual IIA theory. In one dimension we can define a parameter, 

At 

P=~f ( 2 - 7 ) 
which we can think of as a dimensionless inverse temperature, characterizing the behavior of 
the theory. For large (3 » 1 (but still finite as compared with N) the system of DO-branes 
should be well described by a supergravity black hole which is much larger in radius than 
the string length a' 1 / 2 . It is remarkable that since we know how to compute the Bekenstein- 
Hawking entropy of the supergravity black hole, we can predict - assuming the holographic 
correspondence is correct - that in the large j3 limit the precise form of the Yang-Mills entropy 
and free energy will be, 

S = 11.5Af 2 /T 9/5 , / = -4.117V 2 /T 14/5 (2.8) 

where / is the dimensionless free energy, given by F = A 1 / 3 / where F is the usual free energy 
[7]. For decreasing (3 the curvature at the horizon radius becomes larger and the supergravity 
description receives string oscillator a' corrections, about which little is known. 
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For small (3 « 1 the system can best be thought of as a highly excited hot ball of strings 
and branes. Polchinski and Horowitz have argued that the hot ball of strings for [3 « 1 
and the black hole at (3 » 1 are the same object, and the physics at the transition (3^1 
- the 'correspondence point' - is therefore smooth [38]. Witten has argued that the presence 
of a black hole in the dual string theory (geometrically implying a contractible Euclidean 
time circle) indicates the Yang-Mills theory is in a deconfined phase, with thermodynamic 
quantities scaling as 0{N 2 ) and finite expectation value for the amplitude of the Polyakov 
loop < |^Tre^ AdT | > [39]. Conversely the appearance of a confined phase would correspond 
to the absence of a black hole in the dual geometry, or more precisely a non-contractible time 
circle, and results in thermodynamic quantities of order O(l) and a vanishing expectation 
value for the Polyakov loop. Hence at large (3 » 1 we expect the 16 supercharge theory to be 
deconfined as it is indeed dual to a black hole. At small (3 « 1 one can dimensionally reduce 
the theory to a bosonic matrix model, which Monte Carlo simulation shows has energies 
scaling as 0(N 2 ), and hence one expects that the 16 supercharge theory is likely to be 
deconfined for all (3 [33, 34], tallying with smoothness at the 'correspondence point'. 

Hence the key questions that would concern a study of the thermal 16 supercharge Yang- 
Mills theory in the 't Hooft limit would be to confirm the above thermodynamic expectations 
at small and large (3, and to study the transition region to see if there is a phase transition 
in the correspondence region. 

Whilst in this paper we shall not compute with this 16 supercharge theory, we will 
simulate its relative, the 4 supercharge Yang- Mills model in the 't Hooft limit. These two 
theories have similar classical and quantum low energy dynamics and one might expect their 
qualitative thermodynamic properties to be similar. However, strictly one should view the 
calculation in this paper as a warm-up exercise for the 16 supercharge case, and a demon- 
stration that the quantities of interest in the latter theory are likely to be computable using 
similar methods. 

2.2 Infrared behaviour of the 4 supercharge Euclidean theory 

The classical bosonic moduli space of the theory is simply given by setting the scalar and gauge 
field adjoint matrices to be mutually commuting, and constant in Euclidean time. This implies 
they are all diagonal up to gauge transformation. Such a classical moduli space naively leads 
to the concern that the path integral is not well defined due to infra-red divergences. However, 
in both the periodic and thermal cases quantum corrections lift this classical bosonic moduli 
space and hence render the path integral convergent. There are two sources of quantum 
corrections to be considered. These are the off-diagonal elements of the constant modes 
about the time circle, and secondly, the non-constant modes about the Euclidean time circle. 
For a given /3, the loop expansion parameter is determined by the separation of the diagonal 
elements of the constant modes. The infra-red behaviour of the theory is precisely in the 
regime where we take the diagonal elements of the constant modes to be well separated, and 
hence we can approximate the theory by a 1-loop calculation. In appendix B we explicitly 
perform the computation of the effective action for the bosonic classical zero modes in both 
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periodic and thermal cases. We now summarize the results in this appendix, beginning with 
the thermal theory. 



2.2.1 The finite temperature theory 

The thermal case is simpler than the periodic case, and the result we give was previously given 
as a special case in Aharony et al [34]. Taking temperature T, and antiperiodic Euclidean 
time with radius R = 1/T, the thermal boundary conditions imply there are no fermion 
modes that are constant in time, and hence no fermion zero modes. Thinking of the adjoint 
fields as N x N matrices we may expand them as, 



A ab (r) = A«5 ab + -}= £ SA^e% 



2ir 

m=—oo 



. oo 



2vr 

m=— oo 



oc 



= 4= E wWeZ^fr (2.9) 

Z7T 



m=— oo 



where the matrix indices a,b = 1, . . . , N. The classical bosonic moduli are A a , xf. We take the 
perturbations 6A^°\5X^ to have no diagonal terms and since the gauge group is SU (N), the 
field matrices are traceless so the sums 

Ea Ea x i vanish. We define AA ab = R(A a - A b ) 
and Ax" 6 = R(x® — x\), where we note that these are now dimensionless. As we show in the 
appendix, the effective dimensionless coupling constant governing the fluctuations is 

\Ax a 



9^ - r^m ( 2 - 10 ) 



where |Ax af f = ^(Axf) 2 . Hence we may integrate out all the fluctuations with any m 
provided the diagonal component moduli are sufficiently well separated compared to dimen- 
sionless temperature, so, 

\Ax ab \ » /3 3/4 . (2.11) 

As noted in the appendix, this condition arises from the zero modes with m = 0. For small 
circle size, (3 « 1 we can integrate out all 'Kaluza-Klein' modes, ie. bosonic fluctuations 
with m ^ and all the fermionic fluctuations, and are then left with a zero dimensional 
bosonic matrix theory. However, unless the moduli are well separated in the sense above, one 
cannot also integrate out the off-diagonal components of this matrix theory, ie. the m = 
fluctuations. 

As shown in the appendix, when we gauge fix the action and integrate out the fluc- 
tuations at 1-loop we find an action for the bosonic moduli. Supersymmetry is broken by 
the thermal boundary conditions and hence the 1-loop determinants from the bosonic and 
fermionic fluctuations do not cancel each other. We obtain an effective action for the moduli, 

... /cosh|Ax a6 | -cosAA ab \ 

a<b ' 
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This potential takes a simple form, arising from a pair-wise interaction of the moduli and 
hence scales as 0(N 2 ). Gauge invariance implies the potential is periodic in AA ab , the 
moduli A a multiplied by R being angular variables giving the value of the Wilson loop 
p e i§drA _ e if3A a 8 ab _ (X\ag(^e lRAl ) e tRAl , . . . ,e tRAN ^j. The minimum for the potential is 

when the moduli are coincident. Asymptotically, for a large separation |Ax a6 | >> 1 the 
potential goes to zero as 5i_i 00 p ~ — e -\^ ab \ cos /\A ab , and so is attractive towards the coin- 
cident point for -f < AA ab < §. One might worry that all moduli will simply coalesce under 
these attractive pairwise forces. However, as discussed above, for a separation |Ax afe | ~ /3 3 / 4 
the loop approximation breaks down and the theory becomes strongly coupled. Hence we see 
that the classical moduli space is lifted by an attractive potential that drives the infra-red 
dynamics to strong coupling. 

Since the path integral measure over the non-compact zero modes is . dxf we see that 
the integral, giving the partition function, should be convergent in the infra-red due to the 
exponential decay of the potential to zero for large separations. Naively this fast fall off would 
imply that the tails of the eigenvalue distributions of the scalars should die off faster than a 
power law. 

2.2.2 The periodic theory 

The analysis of the periodic case is more subtle. The results we summarize here and give fully 
in the appendix are new, and draw on previous results of Aoki et al and Aharony et al [40, 34] . 
The interesting feature of the periodic case is that in addition to the bosonic classical zero 
modes, there are also fermionic zero modes, where the fermion field matrices are diagonal 
and constant in time. We may then compute a 1-loop effective action for both the boson 
and fermion zero modes by integrating out fluctuations about these. Then integrating over 
the fermion zero modes yields an effective theory for the bosonic zero modes, which has an 
attractive potential that drives the theory to strong coupling. We begin, as above, by writing 
our field as matrices and expanding as, 

oo 

A ab (r) = A«5 ab + J2 SA^e^ 

v m=— oo 
1 oo 

X hah {r) = x?6 ab + -= SXQ***™ 

\J ITS 

v m=— oo 
1 oo 

*a,a b (T) = ejab + ~f= £ ^b^ ^) 
" m=~ oo 

where we have included classical fermionic moduli, We take the fluctuations S^a^ to have 
no diagonal terms and since the gauge group is SU(N), the matrix ^f a is traceless so the sum 

£a vanishes. The dimensionless coupling controlling the integration over the fluctuations 
is again given by g e s as defined above in equation (2.10). 

As shown in the appendix, when we gauge fix the action and integrate out the fluctuations 
at 1-loop we find an action for the bosonic and fermionic moduli. Supersymmetry leads to 
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large cancellations between contributions of the bosonic and fermionic fluctuations, and we 
find that this action is only non-trivial due to the presence of the fermionic moduli, When 
we integrate over these, we then obtain an effective action for the bosonic moduli, 

S!- loop [A a , x 1 }] = - log Yl M aia2 M a2ai . . . M aN - 2aN - 1 M m - iaN (2.14) 

(ai,a 2 ,...,a iv )eP 

where P is the set of permutations of (1, 2, ... , N), and 

oo _. 

M ab = V (2.15) 

m f^ 00 ((m + A^)2 + (Axf)2) 3 

This potential is again periodic in AA ab . It is energetically unfavourable for large separations 
in |Ax" fo |, anc i hence is attractive. As above the classical bosonic moduli space is lifted, 
and the attractive potential drives the theory to strong coupling, with a moduli separation 
|Ax ab | ~ /3 3 / 4 . Since the number of terms in the sum is N\, we expect the 1-loop effective 
action for the bosonic zero modes to have energy 0(N log N). Since the path integral measure 
over these remaining zero modes is rj {dA a Y\ i dxf) we see that the integral, giving the 
partition function, should be convergent in the infra-red. Analogous reasoning to that of 
Krauth and Staudacher [41] in the context of matrix integrals suggests, using naive power 
counting, that the distribution of the bosonic eigenvalues should have tails decaying as a 
power law 1/x 3 . In particular this means that the expectation values of the moments TrX 2p 
for any positive integer p do not exist, even though the partition function itself does. 

As mentioned above, for small enough circle size we may effectively ignore all the non- 
constant modes on the circle and the quantum mechanics should reduce to the 4 supercharge 
matrix integral. This has been studied beyond the 1-loop approximation both numerically 
and analytically. In particular Austing [42] has proven analytically existence of a twisted 
version of the theory and full Monte Carlo simulation was performed by Ambjorn et al [43] 
for large N. 

At this point we should comment that analogous calculations for the 16 supercharge 
theory would yield similar results for the effective potential governing the fluctuations of the 
bosonic zero modes in both the periodic and thermal cases. 

2.3 Phase structure 

In our supersymmetric quantum mechanics with finite N, and hence a finite number of degrees 
of freedom, one cannot have a sharp phase transition, but only smooth cross-over behaviour. 
However, we are interested here in the large N limit of this quantum mechanics. In the infinite 
N limit it is possible to have sharp phase transitions - ie. the cross-over becomes sharper as 
N is increased, leading to non- analytic behaviour for N — > oo. The simplest example of this 
is the Gross- Witten matrix integral which in its large N limit exhibits a 3rd order quantum 
phase transition [44] . Indeed the quenched version of our quantum mechanics theory exhibits 
a discontinuous confinement-deconfinement transition at large N, as discussed in [33, 34]. 
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Hence it is then an interesting question as to whether the supersymmetric theory exhibits 
smooth thermal behaviour or not at large N. As discussed above in section 2.1, holographic 
arguments suggest that the 16 supercharge quantum mechanics is always in a deconfined 
phase, ie. the free energy will scale as 0(N 2 ). Since the infra-red properties of the 16 and 
4 supercharge theories are qualitatively similar one might expect the 4 supercharge theory 
to always be deconfined. Indeed we will see this is borne out in our results, and we see no 
evidence of any sharp phase transition. 

2.4 Ultraviolet behaviour of the theory 

We now consider the UV behaviour of theory, and will show that it is finite. This is extremely 
important for what follows as it allows a naive lattice discretization of the action to recover 
the full supersymmetry of the theory without fine tuning. Let us consider our action (2.1). 
Let us firstly gauge fix the theory so that the gauge field A is constant in Euclidean time, 
d T A = 0. This yields a trivial Jacobian, detd 2 , and allows us then to consider the quantum 
mechanics partition function Z[A] defined as, 



/ 



DAZ[A], Z[A] = J DXi{T)D^{T)D^{T)e~ s ^ ad - Sint (2.16) 



which for fixed A has no gauge dynamics. The unitary matrix A determines the Polyakov loop 
as e lRA , and one performs the matrix integral of Z[A\ over A to compute the full partition 
function Z. The ultraviolet behaviour of the gauged quantum mechanics is determined by 
the gauged fixed quantum mechanics derived from Z[A\. The final matrix integral over A 
does not introduce any new high energy behaviour. 

The quadratic and interaction parts of the quantum mechanics derived from Z[A) for 
fixed A are given by, 

S q uad = Tr j> R dr | X -{D T Xtf + Ma T D T *\ (2.17) 



and, 



," R 

lint = Tr 



[R f 1 A r i2 / \ \ 1 / 2 ~ r i 1 
f dT \-4Nl^M ~{n) *H^'*Jf' ( 2 - 18 ) 



where we have rescaled the fields X{ = (\/N) l l 2 Xi and ^ = (X/N) 1 / 2 ^? to obtain canonical 
kinetic terms. 

Consider now performing perturbation theory in the coupling A, and introducing a high 
energy cut-off A. Clearly quantum mechanics is always super-renormalizable. In general 
this does not preclude divergences, but simply means there can only be a finite number of 
them. Consider the two boson mass renormalization diagrams shown in figure 1. In both 
cases the loop propagators contribute ~ 1/p 2 , and the loop integral J dp/p 2 ~ 0(A _1 ) is 
independent of the cut-off. However for the last diagram in the figure the fermion tadpole 
loop has only a ~ 1/p contribution from the fermion propagator, and superficially the loop 
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Figure 1: Relevant diagrams in super Yang-Mills quantum mechanics perturbation theory. The heavy 
lines are bosonic and the dashed fermionic. The left 2 diagrams are finite mass renormalizations going 
as ~ A -1 . The right diagram is a fermionic tadpole loop which, while naively potentially divergent as 
~ log A, is actually also finite going as ~ A -1 , the divergence not arising due to form of the Yukawa 
gauge interaction. 

integral goes as J A dp/p ~ 0(log A) and hence depends on the cut-off. The only divergences 
in the perturbation theory for this Yang-Mills quantum mechanics come from such fermion 
tadpole loops. 

However, let us consider more carefully the fermion loop in the tadpole. Whilst it superfi- 
cially diverges we will see the gauge and global symmetries actually render the diagram finite. 
We introduce Xf = TiT A Xi and & A = TrT A ty Q where T A are the generators of SU(N) in 
the adjoint representation. The generators obey [T A ,T B ] = if ABC T c with antisymmetric 
structure constants / , and are normalized so that Tr T A T B = 5 AB . Then we may write 
the fermionic interaction term in (2.18) as, 

Cferm ~ f ABC ^a i&a ^X B (2.19) 

The fermion propagator in momentum space is given by, 

1 1 5AB ^ a +0{^) (2.20) 



ia T ^ a D AB 5 6 < a (i5 AB p + f ABC A c ) ip V 

where the subleading terms at large momentum involve the gauge matrix A. The tadpole 
loop for an incoming boson Xf will therefore contribute an integral proportional to, 

fd P (f ABC 6 A B^ + 0(^ =0(1) fdp^O{K^) (2.21) 

where now we see the potentially divergent leading term actually vanishes identically due 
to both the anti-symmetry of the structure constants, / Sab = 0, and also the fact the 
Pauli matrices are traceless so Tr(a*) = 0. Hence we see that actually the fermion tadpole 
loops, and therefore all diagrams in this theory, are cut-off independent provided that the 
momentum regulator preserves the gauge and global symmetry and hence the form of the 
Yukawa interaction term above. This implies that we can regulate our theory with impunity, 
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and providing we have correctly maintained the degrees of freedom given by the quadratic 
action above, and the gauge and global symmetry, we should expect to recover the correct 
theory when we remove the regulator, and in particular all its supersymmetries. This applies 
to the naive lattice regulator we are interested in here. 



/"\ 1/P 

\ 



Figure 2: Relevant diagrams in perturbation theory for Witten's supersymmetric quantum mechanics 
for a quartic superpotential. The left diagram, a mass renormalization, is due to the quartic superpo- 
tcntial term. The right diagram is a fermion tadpole loop due to the cubic term. Both are divergent 
going as ~ log A. 

How generic is this stability against UV radiative corrections? Suppose we wished to 
study the 2 dimensional version of the holographic correspondence we would then have a 
similar 2-d Yang-Mills theory to compute with, and in this case now the mass renormalization 
diagrams in figure 1 above are divergent, since j^dp 2 /p 2 ~ O(logA) as the regulator is 
removed. Unlike the quantum mechanics case there are no internal symmetries preserved 
by a naive momentum regulator that mean the superficial divergence is avoided. Hence a 
naive discretization is not guaranteed to give a continuum limit preserving supersymmetry, 
and one would have to perform a lattice perturbation theory calculation to compute the 
potential counter-terms that would have to be added to counteract this, as for example in 
[45, 46]. Notice that the absence of such counter terms in Yang- Mills quantum mechanics is 
not generic for all quantum mechanics models. For example, Giedt et al [30] have studied the 
case of Witten's supersymmetric quantum mechanics, 

S' = J dr!^x 2 + h\x) 2 + 4>ip + ti'{x)4>^ (2.22) 

Here x is a bosonic field, and ip fermionic, and if we take a superpotential as h(x) = 
mx 2 + 7a; 3 + Xx , then a mass renormalization term is generated as in figure 2 from the 
quartic piece, and the cubic piece gives rise to the fermion tadpole shown in the same figure. 
Naive power counting shows that both superficially diverge as J A dp/p ~ 0(log A), again with 
no internal symmetry preserved by a naive momentum cut-off protecting from this divergence. 
Hence, when discretizing this theory naively one might expect to generate UV sensitive ra- 
diative corrections, and hence break supersymmetry unless the discretization preserves the 
supersymmetry manifestly. Indeed this was confirmed in [47, 30, 29] for the quartic term. So 
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we see that the fact that we are able to discretize our super Yang-Mills quantum mechanics 
is a result not only of its low dimension, but also its interaction structure, in particular, the 
gauge and global symmetry which constrains the boson-fermion interaction vertices which are 
potentially dangerous in quantum mechanics. 

Since we are discussing quantum mechanics we may equivalently use the language of 
operator ordering in order to discuss the UV behaviour of the theory. Starting from a Hamil- 
tonian operator we may construct a lattice path integral discretized in time, by using the 
Hamiltonian to propagate states forward by a small time step in the usual manner. Different 
operator orderings in the Hamiltonian will give different discrete lattice path integrals. If 
there is a physical ordering ambiguity the continuum limit of the discrete path integral will 
yield different physical theories depending on the details of the ordering used. If there is no 
ordering ambiguity, any ordering will yield the same continuum physics. Hence we see that 
in general different naive discretizations of an action will correspond to Hamiltonians with 
different operator orderings, and if there is a physical ambiguity, then different naive lattice 
actions will give different continuum physics. 

We consider again the example of Witten's supersymmetric quantum mechanics. Asso- 
ciated with the continuum interaction term h"(x)^ip in the action is the interaction Hamil- 
tonian, 

H' ferm = Pb h"{x) (2.23) 

with anticommuting fermionic operators b, ¥ so that {¥ , b} = 1. This is the only term in the 
Hamiltonian which is sensitive to operator ordering, as the bosonic interaction term involves 
only x and not its momenta. Indeed we see that changing the operator ordering of the 
fermions changes the purely bosonic interaction terms in the Hamiltonian. For the general 
ordering, 

H'ferm = ((1 " ~ 0) h" (x) + £h"(x) (2.24) 

for some £ we see we have added £h"(x) to the bosonic interaction term. Consider the 
supercharges for this theory, 

Q=(P + ti(xf) S, Qt = P(P-ti(xj) (2.25) 
H'= {&,Q} = P 2 + h'{x) 2 + h"{x) &,b . (2.26) 



and so, 



Hence we see that supersymmetry requires a particular fermion operator ordering. From our 
discussion above, a naive discretization of the action will correspond to a Hamiltonian with 
some operator ordering, and in general it will not be the one required above by supersym- 
metry. Thus naive discretizations will generically lead to continuum theories which are not 
supersymmetric as the continuum bosonic potential will be incorrect. This argument tal- 
lies with our previous considerations of UV behaviour in Witten's quantum mechanics from 
perturbation theory above. 
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Consider now our Yang-Mills quantum mechanics. We again introduce Xf = TiT A Xi 
and ty A = Tr T A ^ a . This fermionic interaction in (2.19) is associated to a Hamiltonian 
interaction term, 

Hferm = f AB °b^m C a Xf (2.27) 
where the fermionic operators St ,b A anticommute as, 

{b^,b*} = 5 AB 5 ita . (2.28) 

Again this is the only term in the action which might be sensitive to operator ordering, as the 
bosonic quadratic interaction term involves only the X{ and not their momenta. However, 

-i,aa f ABC { pA^ $ C } = fABCg^^ = q (2 2Q) 

due to both the antisymmetry of the structure constants and also the traceless property of 
the Pauli matrices. We see that changing the operator ordering in this term introduces no 
additional bosonic terms. From our discussion above, since the form of the Hamiltonian is 
invariant under ordering, naive lattice discretizations will lead to the supersymmetric contin- 
uum physics we are interested in. Thus again we see that while the Yukawa couplings in our 
supersymmetric gauge quantum mechanics superficially could render a naive discretization 
non-supersymmetric in the continuum, the precise structure given by gauge invariance and 
global symmetry ensures that this does not happen. 

We note that since the form of the 16 supercharge action is the same as that of the 4 
supercharge theory, we expect the same argument to ensure that a naive discretization of the 
16 supercharge theory will also regain full supersymmetry in the naive continuum limit. 

3. Two lattice actions 

We now discuss the two actions we have used to simulate the 4 supercharge quantum mechan- 
ics - a naive discretization of the continuum theory and a manifestly supersymmetric lattice 
action. The naive action follows the standard rules for discretizing a gauge theory. We will 
show that the resulting fermion determinant is positive real, which much simplifies simula- 
tion. While the lattice action preserves no supersymmetry, we have argued above that since 
the quantum mechanics is free from UV divergences provided gauge and global symmetries 
are preserved by a UV regulator, all the supersymmetries will be restored in the continuum 
limit - this argument also applies to the 16 supercharge theory of interest for future work. 
The supersymmetric action arises from discretizing a twisted or cohomological formulation of 
the continuum Yang-Mills theory in which nilpotent scalar supercharges can be constructed 
as linear combinations of the original supercharges. In this case we show that half of the 
original supersymmetry of the continuum theory can be preserved in the discrete theory. 
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3.1 Naive lattice action 

In this section we consider a naive discretization of the 4 supercharge quantum mechanics. 
Firstly we rewrite the action using a dimensionless compact coordinate 9, so r = R9, where 
6 has unit period 9 ~ 6 + 1. We create our discrete lattice choosing M lattice sites at 
position 9 n = an, with a the discretization length 1/M. We use R to make our lattice fields 
dimensionless, and denote the variable at lattice site n = 0,1, . . . ,M — 1 with a subscript 
X a , ^ a so that, 

X i {Re n ) = R- 1 X itn , and * a (R6 n ) = R- 3/2 ^ a , n (3.1) 

Since we are taking compact Euclidean time we must identify Xq = Xm, = i^Mj with 
the sign for the fermions giving periodic or antiperiodic (thermal) boundary conditions. 

An important consideration is to avoid fermion doubling on the lattice. This problem 
is encountered if one replaces continuum derivatives by symmetric difference operators. The 
resulting fermion operator possesses additional light states - "doublers" which do not decouple 
as the lattice spacing is sent to zero. Theorems guarantee the appearance of such states in 
any local, translationally invariant and chirally symmetric theory. One simple way to remove 
these unwanted states is to add a so-called Wilson mass term ar A to the lattice action where a 
is the lattice spacing, A is the Laplacian, and r is some non-zero constant. This term lifts the 
doublers to have mass 0(r/a), leaving the physical modes light. For our naive simulations we 
have employed an r = 1 Wilson term, which in quantum mechanics, yields the simple Euler 
discretization prescription ^ a = (^f a +i ~ ^a)/ a - 

We introduce unitary adjoint link fields, U a , to implement the gauge invariance. Gauge 
freedom allows us to fix the link variables on a spanning tree which in our one dimensional 
case means there is one unitary degree of freedom for the whole lattice and we may choose 
how to represent it. For the purpose of simulating the theory a convenient choice is a gauge 
where all the links are equal U = £/j. We may think of U as representing the holonomy of 
the gauge connection about the time circle, so that the Polyakov loop is simply given as, 

e i§Adr = u M_ (3.2) 

Note that the Jacobian introduced by the gauge fixing to set all the links equal is trivial. 
Now we may naively discretize the action as, 

^ (x iA - UXi^U^Y - -^[X iA ,X jA ] 2 + 5^ & ,a (* a ,o " E^a-l^) 

(3.3) 

where again we note the Euler differencing of the fermion kinetic term is free of doublers. 

Let us now show that integrating out the quadratic fermion term gives rise to a positive 
real determinant. As above we again introduce Xf = TiT A Xi and <^ = TrT^^. Taking 
the generators T A to be Hermitian, they obey [T A , T B ] = if ABC T c with the SU(N) structure 



' , — n 



-aa^^^X^] 
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constants f ABC being real. Since Xi are Hermitian the components are real. We then 
define the fermion operator from the above action so that the fermionic part is given by, 
H,aOT AB ^%^d we have, 

aa,AB = § aa ^AB^ _ ^^Tr^C/T^) _ - Q fABC^aa^g (3 4) 

Following a related argument in [43] we see that, 

o*Oi B o* = {0£*Y (3.5) 

and hence any eigenvector v a A a of O having eigenvalue A will be paired with another eigen- 
vector (a 2 ' aa v a Aa)* with eigenvalue A*. Since the determinant of O is the product of its 
eigenvalues, we see it is real and positive. This is very nice as it ensures that in the naive dis- 
cretization we can, in principle, exponentiate the fermion determinant, and use Monte-Carlo 
methods to straightforwardly simulate the resulting action. Notice that we used properties 
of the Weyl representation, and this no longer holds for the 16 supercharge theory. Indeed, 
already in the 16 supercharge matrix integral, the Pfaffian obtained from integrating out the 
fermions is not positive [36]. Thus in the 16 supercharge case, using a naive action one would 
likely have to take the absolute value of the Pfaffian and then use 'reweighting' to simulate 
the phase. How effective this would be would then depend on how important this phase is in 
the physical regime of interest. 

3.2 Manifestly supersymmetric lattice action 

A lattice action which possesses an exact supersymmetry may be derived by dimensional 
reduction of the supersymmetric lattice action for two-dimensional twisted four supercharge 
Yang-Mills theory described in [21, 37]. The resulting action may be written as 

5 = ^E Tr + 2D t<P 3 ) + ^D^W[0 3 ,0] + ^ f [<M +hx^) (3.6) 

where the scalar supercharge arising from the twisting process has the following action on the 
fields 

(3.7) 



QU! 


= ipi 




= ~Dt> 






Q(f> 3 


= K 


Qk 




Qxi 


= B 1 


QBi 


= [4>,x\ 




= V 


Qv 


= [<M 


Q6 


= 



-17- 



It is straightforward to verify that Q 2 acts like a gauge transformation on a generic field / 
which is taken to lie in the adjoint of the gauge group and of the form / = J2iLi ^ f a T a 
where in this section we employ antihermitian generators T a . The Q-exact structure of the 
action then guarantees the discrete action will be supersymmetric. To avoid fermion doubling 
a covariant forward difference operator is utilized whose action on lattice scalar fields is given 
by 

D+(f>{x) = Ui{x)ct){x + i) - <l>(x)Ui(x) (3.8) 

where U\(x) is the usual Wilson gauge link. 1 Notice that several of these fields, such as U\(x), 
carry a vector index. In the continuum this is redundant for a one dimensional theory but in 
this lattice construction the index plays an important role, indicating that such fields live on 
the links of the lattice and transform under gauge transformations as 

f l {x)^G{x)f l {x)G\x + l) (3.9) 

where G(x) = e^(x). The corresponding transformation for scalar or site fields is 

f(x) - G{x)f{x)G\x) (3.10) 

These transformation reduce to the usual continuum ones in the naive continuum limit. 
Notice also that the definition of the forward difference operator when acting on a site field 
automatically produces a vector or link field with the correct lattice gauge transformation 
properties. 

Carrying out the Q-variation and integrating out the auxiliary field B\ we find the fol- 
lowing lattice action 

S = Yl Tr f " Yl xiD ~ D+xi ~ X ^ ~ X i D tp i - P ijD I X i + 5 yk J (3.11) 

x \ i=l i>j J 

where we have written (f> = X 1 + iX 2 , <p = X 1 — iX 2 , cf> 3 = X 3 and the covariant backward 
derivative , which is adjoint to D^, acts on link fields in the following way 

D^tpi = %b!( x )U\x) - U\x)^i{x - i) (3.12) 

We have also relabeled the fermions according to 

\\ = Vi (3.13) 
A? = xi 

(? = K 



1 actually a complexified version of it since the construction requires the lattice fields be taken complex 
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Notice that the bosonic action is real positive definite on account of the antihermitian basis 
for the fields. The Yukawa interactions take the form 



s YK = p 2t [*V] - *P 2t [*V] - P lf [*V] - V f [*V] 

+ \»[X\ X\] - *A^[X 2 , A}] - \?[X\ A 2 ] - iA? [X 2 , A 2 ] 
+ {A^A 2 ]-p 2 t[X 3 y]+h.c} 



(3.14) 



While the bosonic action arising in this discretization is rather similar to the naive action 
described in the previous section, differing mainly in the definition of the lattice derivative, 
it is clear that the fermionic action is quite different - the fermions can be assembled into 
a four component object which in the continuum is just the usual Majorana fermion of four 
dimensional N = 1 super Yang-Mills. Thus the result of integrating out the fermions in this 
twisted formulation is to produce a Pfaffian rather than the determinant encountered with 
the naive discretization. Of course, in the continuum, one can find a change of variables 
which allows the Pfaffian to be rewritten as a simple determinant but this is no longer true in 
the lattice construction. Indeed the complex two component spinor encountered in the naive 
formulation takes the form 



Since the different component fields carry different gauge transformation properties this spinor 
does not transform simply under lattice gauge transformations and cannot be used to con- 
struct a gauge invariant lattice action. 

A more difficult question relates to the complex nature of the fields in the supersymmetric 
lattice construction. To target the correct continuum theory requires choosing the path 
integral along a contour such that the imaginary parts of all fields are zero. The question 
then arises as to whether the supersymmetric Ward identities are satisfied along this path. 
The latter are Q-exact and hence by standard arguments can be evaluated exactly in the 
semiclassical limit in which the fields are expanded about the classical solution [26] . It is not 
hard to show that in this limit the real and imaginary parts of the fields decouple and the 
computation can be consistently truncated to the real line. These theoretical arguments are 
strengthened by the results of Monte Carlo calculations which support the existence of an 
exact Q-supersymmetry at the quantum level [37, 48]. 

Finally notice that this lattice action actually possesses a global symmetry of the form 




(3.15) 
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This symmetry combined with the original Q-supersymmetry leads to the lattice action pos- 
sessing a second exact supersymmetry. The existence of this second supersymmetry is con- 
sistent with the orbifold construction described in [13]. 



4. Results 

In this section we discuss the results of lattice simulation of our two implementations. We 
begin by considering the quenched theory, move onto the periodic theory and end the section 
by discussing the thermal theory. The results for the quenched theory agree precisely for the 
naive and supersymmetric actions since they differ only in their treatment of the fermions 
(notice though that no gauge fixing is done in the supersymmetric implementation). As we 
discuss, the naive and supersymmetric implementations both give compatible results for the 
periodic theory. In practice the naive implementation is computationally easier since the 
corresponding fermion operator is half that of the supersymmetric formulation and, as we 
have shown, is real and positive. Hence the bulk of the periodic results were generated with 
the naive action. In the case of the thermal system the situation is more interesting - for 
coarse lattice spacings and a small number of colors we observed that both codes suffered from 
strong lattice artifacts most visible in the distribution of scalar eigenvalues which developed 
long tails out to large eigenvalue - a situation quite different from the complementary runs 
with periodic boundary conditions. We observe that this stems from very large fluctuations of 
the scalars in the classical moduli space - for the periodic system the corresponding bosonic 
zero mode is strongly suppressed by its fermionic superpartner but in the thermal case the 
would be fermionic zero mode is lifted by the antiperiodic boundary conditions and is less 
effective at finite (3 at suppressing these zero mode effects. This problem appears to be worse 
for the supersymmetric action and hence we have again derived the bulk of our thermal results 
from the naive action runs. In the latter case these effects appeared negligible for N > 5 and 
M > 5. 

We have examined the following expectation values. We normalize the expectation value 
of the fermionic action as, 

&f = 2M(N 2 - 1) < ^ fcrmionic > ( 4 " 1 ) 

and a simple scaling argument gives the Schwinger-Dyson equation Sp = 1 which we then 
use to check the equilibration of our runs. We also compute the bosonic action as, 

Sb = ^2 (< ^bosonic > ~S ze ro) (4-2) 

where we have subtracted the extensive zero point energy contribution, which is given ex- 
plicitly as |M(iV 2 — 1) for both actions as can easily be seen in the weakly coupled high 
temperature limit. For the supersymmetric lattice action we have Sb = for all (3 since it is 
related to an index [21]. For the naive periodic theory we expect Sb to approach zero as the 
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number of lattice points M — > oo. We have normalized the bosonic action by N 2 to ensure 
that Sb should tend to a constant at large N in a deconfined phase. 

This form of the bosonic action has an additional interpretation in the thermal and 
quenched theories as yielding a measurement of the dimensionless mean energy of the system, 
< E >. This dimensionless energy is given by the usual relation < E >= — 9 q^ Z so that, 

S B = -| < E > . (4.3) 

Thus the vanishing of (the subtracted) Sb in the periodic theory corresponds to the usual 
requirement that the supersymmetric theory have vanishing vacuum energy at zero temper- 
ature. 

We are also interested in the behaviour of the Polyakov loop variable, and we compute 
the expectation value of its modulus, 

P= -L < \Tre^ AdT \ > (4.4) 

and also its corresponding susceptibility 

dP = ~< \ ^Tre* f Adr ^j 2 | > - < \Tve^ AdT \ > 2 . (4.5) 

We note the inclusion of 1/iV in these definitions to ensure that at large TV" in a deconfined 
phase, P, dP should tend to a constant. 

We have also computed the distribution of eigenvalues of the scalar fields, made dimen- 
sionless as RXi(r), and averaged over the lattice. We have seen no evidence of broken 50(3) 
symmetry 2 and so assume the distribution for i = 1,2,3 is the same, and denote it as x(//). 
We observe this appears to have a well-defined large TV-limit with a width which is controlled 
by (3. From our earlier discussion of the 1-loop eigenvalue potentials we deduced that well sep- 
arated dimensionless eigenvalues were attracted together until they entered a strongly coupled 
infra-red regime when their separation goes as ~ /3 1//4 . This therefore is the expected scale 
characterizing the parametric dependence of the width of the distribution x(/x), and indeed 
this is borne out by our results later. For SU(N) the distribution has N peaks localized close 
to the origin, and as we have discussed we expect the distribution to fall off most slowly for 
the periodic theory, going as x(/jl) ~ H -3 for large fi. This tail appears to be universal for 
the periodic boundary conditions depending on (5 but not N [48]. The tail means that the 
standard deviation is in principle ill-defined as an observable, at least for periodic boundary 
conditions. Therefore in order that we may characterize the quenched, periodic and thermal 
theory scalar distributions we define the width by the observable, 

W = /T 1/4 J dfj, \n\ x(jj) (4.6) 

which is well defined even for the expected tail behaviour of the periodic theory. 

2 We note that at finite N such a spontaneous symmetry breaking is not possible, although it might be in 
the large iV limit. 
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4.1 The quenched approximation 

We begin by discussing the quenched theory, where the fermions are simply ignored. The 
lattice simulation of this bosonic theory is then a very tractable problem, and one can easily 
work at large ./V and establish the 't Hooft scaling. Such gauged quantum mechanics with 
adjoint scalar fields was originally studied in [33, 34] where evidence was given that the theory 
undergoes a large N confinement /deconfinement transition. We now review this behaviour 
as it indicates where the interesting dynamics in the unquenched theory is likely to occur, 
and how we might see this in the available observables. 

The first observation we may make from the quenched theory is that the continuum limit 
is very easy to obtain. In the appendix D we show the bosonic action and Polyakov loop for 
different numbers of lattice points M = 5,8, 12. We see that the results we obtain for all 
these quantities on a lattice with M = 5 are very close (within one percent) to those with 
M = 8, 12. Indeed we find this to be true for all the observables we measure in the quenched, 
periodic and thermal theories (with the exception of the bosonic action for the naive periodic 
theory which we discuss later). Hence, in the main text, results are calculated using M = 5 
lattice points unless otherwise stated. Certainly, in the case of the unquenched theory, the 
statistical errors in most of our measurements are larger than the systematic discretization 
error incurred by using the small lattice for the range of (3 studied. 

In figure 3 we plot the expectation value of the bosonic action Sb, the expectation value 
of the modulus of the Polyakov loop, P, and its corresponding susceptibility dP for various 
numbers of colors N up to 16. Indeed it would be easy to compute at larger N but this is not 
our objective here. As observed earlier in [33, 34] we see that there appears to be a sharp, 
probably first order large N phase transition in behaviour at f3 = (3 C ~ 0.85, signaled by the 
bosonic action Sb, the Polyakov loop P, its susceptibility dP and also the scalar eigenvalue 
distribution width W. In particular, the Polyakov loop variable P remains finite for (3 < f3 c , 
but appears to be consistent with zero at large N for (3 > (3 C (see [33, 34] for data with N up 
to 30). Furthermore, the mean energy < E > appears to suffer a discontinuity around the 
critical point. 

Notice also that the linear regime observed in Sb at large (3 is consistent with a non-zero 
vacuum energy - the latter being given by the slope of Sb = —(3/3 < E > with respect to (3 at 
large (3. This is, of course, consistent with what one would expect for a non-supersymmetric 
theory. Conversely the observed /3-independence of Sb at high temperature (small (3) is 
consistent with the expected classical thermal behavior < E >~ N 2 T for a theory possessing 
iV 2 deconfined gluons. 

4.2 Full supersymmetric theory: Periodic 

We now focus on the theory with dynamical fermions having periodic boundary conditions 
- ie. the theory in finite volume and zero temperature. Our task is to compare the results 
obtained from our two independent implementations - the naive and supersymmetric - and 
show they are consistent, and produce the correct supersymmetric continuum physics. 



- 22 - 




Figure 3: Plots of action Sb, Polyakov loop P, its susceptibility dP and the scalar width W for the 
quenched theory for varying N, and 5 lattice points. 



Figure 4 shows a comparison of the Polyakov loop observable P and its susceptibility 
dP for the two lattice implementations in the case of N = 3 and M = 5 lattice points. We 
see the results are in in rather close agreement, in fact within statistical error, for P, dP. 
The close agreement is a good check on both implementations as it indicates that lattice 
spacing effects appear to be small and supports the claim that the M = 5 runs with naive 
action yield results which are already decent approximations to the continuum. The figure 
also shows a detailed comparison of the scalar eigenvalue distributions for the naive and 
supersymmetric actions and again reasonable agreement is seen. We study the continuum 
behaviour of the naive periodic theory in appendix D, and confirm that the observables P and 
dP are close to their continuum values by simulating at M = 8 and 12 lattice points. Given 
constraints on resources we have not been able to compute for larger N in the supersymmetric 
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Figure 4: Top: Plot comparing scalar field distributions for the 2 implementations with periodic 
fermions (bars for naive, dots for supersymmetric; left (3 3 = 0.1, right (3 3 — 2.0). Bottom: the 
Polyakov loop P and its susceptibility dP for N — 3 and 5 lattice points. 

implementation and therefore cannot check the agreement there. 

We show the action Sb for the supersymmetric implementation in figure 6 for N = 5 and 
5 lattice points and see that it is consistent with zero, as we expect for the vacuum energy of 
a system with exact supersymmetry. 

Figure 5 shows the bosonic action Sb for the naive theory at three representative values of 
(3 s with N = 3,5 as a function of increasing numbers of lattice points. While by construction 
this vanishes in the supersymmetric theory, in the naive theory we only expect it to vanish 
in the limit M — ► oo. We see, both for N = 3 and 5 and fixed /3 3 , that the numerical results 
are indeed consistent with Sb decreasing to zero as M increases. The time taken to compute 
larger lattice sizes with reasonable statistical errors has prohibited extending these plots to 
higher numbers of lattice points. While the approach to the continuum value of the action 
is rather slow, appearing to go as ~ 1/M P with p ~ j ~ g> even fo r the modest number of 
lattice points M = 5 its actual value is already rather small compared to the quenched or 




Figure 5: Plots of the bosonic action Sb for the naive discretization with periodic fermion boundary 
condition for N = 3 (left) and N = 5 (right) for increasing lattice size M = 5, 8, 12 for 3 values of 
/3 3 = 0.5,1,2. 
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Figure 6: Plot of the bosonic action Sb for the supersymmetric discretization with periodic fermion 
boundary condition for N = 5 with M — 5 lattice points. We see it is consistent with zero as we 
expect. 



thermal theory (eg. see figures 3 and 8) for these values of (5. 

In figure 7 we plot the Polyakov loop P, its susceptibility dP and the scalar eigenvalue 
distribution width W for N = 3 compared with those of the quenched theory. Interest- 
ingly, while the action and scalar distribution width are different to the quenched theory, the 
Polyakov loop behaves rather similarly. In figure 8 we also see the same is true for N = 5. 
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Figure 7: Plots of the Polaykov loop P, its susceptibility dP and the scalar width W for the periodic 
theory for N = 3 with 5 lattice points, also compared to the quenched theory. 



We have currently not simulated the periodic theory for N > 5 and therefore cannot confirm 
this correspondence with the periodic theory occurs in the large N limit. If it does, it leads 
to the interesting conclusion that while the bosonic action Sb is a constant in the periodic 
continuum theory independent of volume and coupling, the behaviour of the Polyakov loop 
may nevertheless not be smooth in the large iV limit. We note that the scalar eigenvalue 
width is broader for the periodic theory than for the quenched. This is consistent with the 
1-loop calculations of the potentials on the classical moduli space where we expect the peri- 
odic theory will have power law tails in the eigenvalue distributions, whereas the quenched 
will not. 

Let us summarize these results. We have compared the supersymmetric and naive lattice 
actions for dynamical fermions with periodic boundary conditions. The results are consistent, 
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and both actions yield the expected behaviour with the supersymmetric version giving a 
vanishing bosonic action and the naive version giving a continuum limit consistent with 
vanishing bosonic action. Hence our earlier analytic claims that a naive discretization of the 
action will give the correct supersymmetric continuum physics appears to be borne out in 
practice. 

4.3 Full supersymmetric theory: Finite temperature results 

We now turn to the most interesting part of our results: the supersymmetric theory at finite 
temperature. If we were studying the 16 supercharge theory these lattice simulations would 
be dual to a computation of the thermal properties of N DO-branes, and should reproduce 
the thermodynamics of black holes described earlier at low temperature. In the case of the 
4 supercharge model studied in this paper no such correspondence exists. Nevertheless, we 
might expect on the basis of the one loop calculations described earlier that this model lies in 
a similar universality class and so understanding how to extract continuum results from this 
model should stand us in good stead for a future simulation of the thermal 16 supercharge 
theory. 

We have observed that the thermal theory exhibits large lattice artifacts for small values 
of N and the number of lattice points M . These lead to an apparent instability in the scalar 
eigenvalues. The problems seem most acute with the supersymmetric action and so we have 
concentrated on using the naive implementation for the bulk of our thermal runs. In this 
case to avoid these strong artifacts we require N > 5, M > 5 for the range of (3 we are 
studying, /3 3 < 10.0, and then a good lattice continuum limit is seen. We discuss these effects 
in more detail in the appendix D, and it is an interesting direction for future research to 
better understand how to more stably implement the supersymmetric lattice action at finite 
temperature. 

We begin in figure 8 by comparing this thermal theory with the previously discussed 
quenched and periodic theories for N = 5 and M = 5 lattice points. In appendix D we show 
these quantities for increasing lattice sizes, showing that already these 5 lattice point data 
capture the continuum reasonably accurately. 

Our expectation for this theory is that at small (3 (high temperature) the thermal theory 
will behave like the quenched theory since the fermions are lifted out of the dynamics by their 
thermal mass. On the other hand in the absence of spontaneous supersymmetry breaking, in 
the large (3 limit (low temperature) we expect the energy of the theory to go to zero and the 
behaviour to coincide with the periodic theory. 

We clearly see from the data that at small (3 the action Sb, Polyakov loop P and scalar 
eigenvalue distribution width W do coincide with the quenched behaviour as expected. Con- 
versely, at large (3 we see the same observables depart from the corresponding quenched 
quantities and approach the periodic results. The variation of the bosonic action, Sb, is 
flatter and apparently smoother than that of the quenched theory and appears to have a 
small linear slope for large (3, which would imply a non-zero energy and hence supersymme- 
try breaking if it were to continue to large values of (3. However, given we have data only 
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Figure 8: Plots of the bosonic action Sb for the quenched and thermal theories, and the Polaykov 
loop P and scalar width W for the quenched, periodic and thermal theories for N — 5 and 5 lattice 
points. These quantities are expected to be close to their continuum values. 



to inverse temperatures (3 ~ 2.2 it is difficult to say whether we are really seeing asymptotic 
behavior - any sublinear behavior would imply supersymmetry is restored at zero tempera- 
ture. And the apparent asymptotics may also be influenced by discretization effects. Indeed, 
to address the possibility of supersymmetry breaking one should first extrapolate the data to 
zero lattice spacing at fixed (3 and then examine the (5 dependence of the extrapolated curve. 
This extrapolation is beyond our current resources. Nevertheless, it should be noted that the 
asymptotic behavior of this quantity appears to be rather different than that expected for the 
16 supercharge theory where we see from equation (2.8) that holography predicts Sb ~ /3 -9 ^ 5 . 

Independent of the final conclusion concerning supersymmetry breaking there is little 
indication of a discontinuity in the mean energy as a function of (5 as N increases and our 
numerical results are consistent with the existence of just a single phase for all /?. In figure 9 
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we show the bosonic action Sb, Polyakov loop variable P, its susceptibility dP and the scalar 
eigenvalue width W for the thermal theory for N = 5, 8, 12. We see the results confirm the 't 
Hooft scaling - recall we expect the observables Sb, P, W to tend to a constant at fixed (3 in 
the large ./V limit, as we see confirmed in the data. Whilst in the quenched theory there is a 
large N transition, with P vanishing at large N for (3 > (3 C , and the susceptibility diverging 
at /3 = (3 C , we see no such behaviour here for the thermal theory in the range of (3 studied. 
Instead P appears smooth and the susceptibility dP reduces for increasing N over most of 
the range of (3. Together with the apparently smooth bosonic action Sb and scalar width 
W we conclude that the thermal theory is always in one, presumably deconfined, phase with 
non-vanishing P and with bosonic energy 0(N 2 ). 




0.5 1 1.5 2 0.5 1 1.5 2 




0.5 1 1.5 2 0.5 1 1.5 



Figure 9: Plots of the bosonic action Sb, Polaykov loop P, its susceptability dP and the scalar field 
width W for the thermal theory for increasing N — 5, 8, 12 with 5 lattice points. These quantities are 
expected to be close to their continuum values and show no indication of large N phase transitions. 
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5. Discussion 



This paper is devoted to a study of four supercharge Yang Mills quantum mechanics at large 
N. This work is motivated by the idea that a related model - namely the sixteen supercharge 
theory, should be dual to IIA string theory at least at sufficiently low energies that it can 
be approximated by a supergravity theory In this regime the Yang-Mills model is strongly 
coupled and hence we have developed lattice discretizations of the model which allow for 
Monte Carlo simulation. Two such discretizations have been studied - a so-called naive 
action in which supersymmetry is broken classically by terms of the order the lattice spacing 
and a manifestly supersymmetric action which arises from a discretization of a twisted form 
of the continuum theory 

We have analytically argued that the quantum mechanics is independent of a UV mo- 
mentum regulator provided it preserves the gauge and global symmetries and hence our naive 
discretization should give the correct continuum physics. The mean energy and the Polyakov 
loop computed with these two discretizations agree well for the case when supersymmetry 
preserving periodic boundary conditions are used for the fermions, and the mean energy of 
the naive discretization is indeed consistent with vanishing in the continuum as required by 
supersymmetry. 

The theory has a large classical bosonic moduli space. We have analytically computed 
the 1-loop effective action for the bosonic moduli in this theory - both for periodic and 
antiperiodic (thermal) fermion boundary conditions and find that the classical moduli space 
is lifted and the theory does not suffer from infrared divergences. These analytical results are 
in agreement with our numerical work which shows that the eigenvalues of the scalar fields are 
localized close to the origin in field space where the low energy behaviour is strongly coupled. 

We have studied the thermodynamic behaviour of the thermal theory for N < 12 and ob- 
serve the expected 't Hooft scaling of thermodynamic quantities. In contrast to the quenched 
theory, the behaviour of the thermal theory with dynamical fermions appears smooth and we 
find no evidence of a large N phase transition, the theory always appearing deconfined. This 
is similar to holographic expectations for the 16 supercharge theory, although we find that 
the temperature dependence of the free energy in the 4 supercharge model is rather different 
from that expected for 16 supercharges. 

It would be very interesting to extend these calculations to the case of 16 supercharges. 
This is computationally more challenging with the main problem being how effectively the 
Monte Carlo procedures can handle the phase of the Pfaffian arising after integration over 
the fermions in that case. This problem is currently being studied. If calculations with this 
Pfaffian prove possible it would be extremely interesting as then the lattice model could be 
used to further test the duality between gauge theories and gravity and perhaps learn more 
about the nature of the gravitational theory in regimes of high temperature where stringy 
corrections are not small. 
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A. Holographic dual of 16 supercharge quantum mechanics 

Following Itzhaki et al [7], we consider the "decoupling" limit of N coincident DO-branes. 
We take N large with Ng s fixed, where g s is the string coupling. The decoupling limit is 
then defined by considering excitations of these DO-branes with fixed energy while sending 
the string length scale to zero so a' ^ 0. In this limit the degrees of freedom of the system 
split up into those localized near the branes - the 'near horizon' excitations - and those living 
far from the brane which we are not interested in here. 

Depending on Ng s there are two perturbative descriptions of the degrees of freedom 
living near the branes. For Ng s << 1 the DO-branes decouple from the ambient 10-d gravity 
and the degrees of freedom are well described simply by the worldvolume theory of the DO- 
branes whose degrees of freedom are the open strings ending on the branes, with dynamics 
governed by the 16 supercharge SU(N) Yang-Mills quantum mechanics. Since we consider 
fixed energy excitations and a' = Z 2 . — > 0, the action is just the conventional two derivative 
one, the higher derivative a' corrections being irrelevant. The gauge coupling is then found 
to be gyM = 9so'~ s ^ 2 /(2-7r) 2 and in this Ng s « 1 regime it is small. 

For Ng s >> 1 the DO-branes couple strongly to gravity and the appropriate perturbative 
description is given by the target spacetime supergravity solution for N DO-branes. The 
string frame metric is 

ds} IA = -^1 + ^/1 Qdr 2 +r 2 ^ 2 8) ), 

/(r) = l-(^) 7 , h(r) = l+(^) (A.l) 

and the dilaton and Ramond-Ramond (RR) gauge field are given as, 

e* = g s hi, A ^ = -l(l- 1 )- (A-2) 
The solution describes a black hole with temperature T and Bekenstein-Hawking entropy S, 

using the conventions of [7] with Os the area of a unit 8-sphere. Each brane carries RR 
charge, and hence the number of DO-branes is computed from the charge of the black hole by 
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N = J ss *F with F the RR field strength, and one finds, 



d7 / r 7\l/2 

N ^k( 1 + w) ■ < A ' 4 > 



Hence in the limit Ng s is large the characteristic curvature radius of the solution, estimated 
by R, is much larger than the string length and therefore the string worldsheet theory is 
weakly coupled, and supergravity is a good approximation. Geometrically this solution has 
two regions. One is asymptotically flat 10-d spacetime, the other is the geometry describing 
the 'near-horizon' region of the branes. Finite energy excitations in these two regions are 
separated from each other by a potential barrier. To take the decoupling limit and focus 
on the excitations only in the 'near-horizon' region we must fix our physical energy scales of 
interest and take a' — > 0. We identify the energy scales we wish to fix as U and Uq, where 

U=-„ U =^ (A.5) 
a' a' 

and so U > Uq. Then taking a' — > keeping U, Uq fixed we see that Ng s — ► ^ and the 
near- horizon geometry becomes, 

7 



f(U) = 1 - ^ (A.6) 
with A = Ngy M . In this decoupling limit the entropy becomes, 

and the temperature is given by, 

r/Al/3 = 16 7n^( c/ ° /Al/3 ) 5/2 ' (A8) 



The key statement of Maldacena's duality is that both these perturbative desciptions - 
the Yang- Mills for Ng s « 1 and the stringy black hole for Ng s » 1 - are not limited to the 
regime of Ng s where they are perturbatively good. In particular the Yang-Mills description 
remains well defined for all Ng s . The string black hole description also persists for finite Ng s 
away from Ng s » 1, and for Ng s finite one must take into account stringy a' corrections to 
the description of the black hole in supergravity. 

To summarize the correspondence, the claim is that Yang-Mills with coupling gyu at 
finite temperature T taken in the large iV limit is dual to string theory with target space 
given in the large Ng s limit by the near horizon geometry of N DO-branes as above, given by 
fixing our thermal energy scale Uq and taking a' — > 0. The dimensionless effective coupling at 
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finite temperature in the Yang-Mills is given by (3 3 = X/T 3 . The equations (A. 7) and (A. 8) 
then relate the Yang-Mills quantities to the string theory. 

Let us further explore the corrections to the supergravity description of the string theory 
above. It is crucial that all curvatures and the dilaton are small in order that the above 
supergravity solution is valid. The curvature radius p at energy scale U is characterized by 
the radius of the sphere in the above geometry, so that in string units, 



P ( A 



a 



ll/2 \ JJ3 



and the dilaton at the radius U is, 



-(-) 
N\U 3 



1/4 

(A.9) 



7/4 

(A.10) 



We also require the temperature T to be large enough that Euclidean winding modes are not 
present at the horizon. 

Hence we see that provided X/U 3 , X/Uq and N are large the supergravity solution above 
is a good approximation. In particular the dilaton condition shows we must take the 't Hooft 
limit, first taking N to infinity with X/U 3 , X/Uq fixed, and then take these large. 

If we take the black hole to have a high energy /temperature with X/Uq ~ 1 we reach 
the Horowitz-Polchinski correspondence region where the IIA supergravity breaks down even 
at the horizon due to a' curvature corrections becoming important. One requires full IIA 
string theory to describe the horizon region. Conversely for a black hole with ultra low 
energy/temperature, so Uq/X ~ iV -4 ' 7 , which is outside the 't Hooft scaling limit, the dilaton 
becomes large near the horizon and the string theory becomes strongly coupled there. It 
may then be resolved by lifting to M-theory where again one finds 11-d supergravity is valid. 
However, we emphasize that provided we remain within the 't Hooft scaling regime, we cannot 
access such low temperatures. 



B. Potentials on classical moduli space 

We now derive the effective 1-loop actions for the classical zero modes. We compute first 
the potential for the quenched theory, and derive the effective coupling controlling the 1-loop 
integration. Then we proceed to the thermal theory, and finally the periodic theory. 

B.l Quenched case 

The results in this subsection are a special case studied in Aharony et al [34] . We expand our 
bosonic fields in fluctuations about the classical zero-modes, 

A ab (r) = A a 5 ab + 5A ab (r) 
Xi, ab (r) = x?6 ab + 8X i>ab (T) (B.l) 
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where the the classical bosonic moduli are A a ,xf. We take the fluctuations SA,6Xi to have 
no constant component in their diagonal terms. We will shortly expand these fluctuations in 
harmonics so, 



^ oo 



'2tt 

m=— oo 



1 °° 



and thus 5A^,<5X^ fe vanish when a = b. Since -A^r), X itab (r) are Hermitian fields, we 

have (6A^)* = SA^ m ^ and (SX^)* = $X^™\ Since the gauge group is SU(N), the sums 
J2 a A a , Y, a xf vanish. We define the dimensionless AA ab = R{A a -A b ) and Axf = R(xf-x%), 
and the following quantities, 

5®/j,,ab = (5A ab ,5Xi tab ) 

Df = (ltd T + iAA ab , iAxf ) (B.3) 
We may now write the action to quadratic order in the perturbations, obtaining, 

squad = ^ | dr £ S% >ab Mp^ ab (B.4) 

a<b 

where, 



Mf b 
n» 2 



\ (DfDf - 5, v {Dff) . (B.5) 



We must take care with the action since we have not fixed a gauge. The gauge transfor- 
mation generates a bosonic fluctuation 8<&^ ab = D°^\ ab for arbitrary functions X ab of r and 
correspondingly we see, 

Mf u D? = Va,6 (B.6) 

We choose to use the gauge freedom to fix the gauge 5A ab {r) = 0, and hence the gauge field 
is diagonal and constant in r. This gauge fixing introduces a Jacobian factor of \\ a b det D® b 
from the path integral measure. In the case of the constant modes on the circle, this is just 
the familiar Vandermonde determinant. Now the quadratic action becomes, 

squad = JLJ dT ^2 (sxt^MgSX^) - J> I det Df\ (B.7) 

a<b a,b 

We may now expand in Fourier modes to obtain, 

- j 00 00 

squad = ^Y.H {sX^M^SX^y) - £ £ InpW (B.8) 

a<b m=—oo a<b m=— oo 
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where D { ™ )ab = (2vrim + iAA ab , iAxf), and 

M^ ab = l - (^ D ^b D bn)ab _ s ^ D (m)ab^ ( B .9) 

For a given mode m, and colour indices ab, the eigenvalues of the 3x3 matrix M^ ah are 

—{D < ^ ab ) 2 once, and — ^ a {D^ ab ) 2 repeated twice. Now our quadratic action is positive 
after our gauge fixing, we may integrate out the fluctuations yielding, up to terms independent 
of the moduli, 

oo 

gquad _ 2 'y ^ 'y * J n |^)(m)ab|2 

a<b m=— oo 

oo 2 

= 2^ Yl ln((2vrm + A^ ab ) 2 + |Aa; afe | 2 ) (B.IO) 
a<b m=— oo 

where |Ax a6 | 2 = ^2i(Axf b ) 2 . We note that the constant m = modes have an enhanced 
50(4) target spacetime global symmetry. This is explicitly broken by our gauge choice, but 
the result should have this full invariance, and we indeed see that this is the case. We also 
note that the effective action takes the form of a pair-wise interaction between the moduli 
with colour labels a and b. 

The infinite sum over Fourier modes must be regulated in the quenched case. A Pauli- 
Villars regulator can be utilized, giving a regularized 1-loop action, 

/(2vrm + AA ab ) 2 + |Ax ab | 2 ^ 2 



a<bm=— oo \ \ / / 

= 2 Yj 1u ( cosll \Ax ab \ - cosA^ ab ) (B.ll) 



a<b 

where in the second line we have suppressed the trivial divergence going as R£l. 

When was this 1-loop integration valid? As with any dimensional reduction on a circle, 
the strongest coupled modes are the ones that are constant on the circle. We should compare 
our quadratic action (B.8) for these constant modes with the interaction terms, 

= A^ RTl ( 5**(o)ff*(p)ff**(o)(y*(p) - S^S^S^S^A (B.12) 
2 A V / 

In order to obtain a canonical kinetic term in (B.8) we should rescale the constant fluctuation 
fields, 

54 <o> /^5 w <o> (R13) 



Our interaction term now takes the schematic form, 

A-R 3 1 , .--ffVKA 1 



* mt ~ ^Ia73 ( ^ (0))4 = ^iv3 ( ^ ( ° ))4 (B - 14) 
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where we have suppressed the obvious index and matrix structure. The factors of N are 
associated with the 't Hooft limit, leaving the effective coupling for the constant modes, 



/3 3 

9eff = o (B.15) 

U ((AA ab ) 2 + \Ax ab \ 2 f 

Since RA a should be thought of as an angular variable the effective coupling is really char- 
acterized by just the non-compact bosonic moduli, 

/3 3 

9eff ~ (B.16) 
and hence we expect the 1-loop approximation is good provided, 

|Ax afc | » /3 3 / 4 . (B.17) 

B.2 Thermal case 

Again this case is straightforward, and closely follows the quenched calculation above. The 
results in this section are a special case studied in Aharony et al [34]. We expand our fields 
as before, now including the fermions with the appropriate boundary conditions, so, 

oo 



A ab (r) = A«5 ab + J= £ SA^e^ 

\/ ZTT 

v m=— oo 
1 oo 



/2tt 

v m=— oo 
1 oo 

= 4= E W^e&^fr (B.18) 

v m=— oo 

We see that due to the antiperiodicity the fermions can have no constant mode on the circle. 
As before we take 8A$ , 8Xf^ and now also S^S^ b to vanish when a = b. We gauge fix the 
action in the same manner, and now arrive at a similar action to the quenched one above in 
equation (B.8), but with a fermionic piece too, 



AT 00 
gquad __ ^ ^ 



E {sx^m^&xQ) + e W m+I)a6 ^) - E ln \ D i m)ab 



a<b a,b / a<b 



(B.19) 

Now performing the fermionic integral yields a determinant equal to rim=-oo EL b (D^ + ^ ab ) 2 ■ 
The bosonic fluctuation integral yields the same result as for the quenched theory. Putting 
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this together we obtain the effective action for the bosonic zero modes, 



S quad = 2 Y^ ^ ln\D^ ab \ 2 ln \ D 

a<b m=— oo a,b m=—oo 



(m+i)a6,2 
a 



_ 2 y- y> ( (27rm + A^ ab ) 2 + |A3; Qfe | 2 \ 2 
~ 2^ 2^ n I (2vrm + vr + AA ab ) 2 + |Ax afe | 2 / 

a<bm=— oo xv / i i / 

_ 2 V^i /cosh |Ax a6 | - cos AA ab \ fB20l 
^ Vcosh |Ax a6 | + cos AA ai> / 

where no regulator is required to evaluate the infinite sum in the last line. 

B.3 Periodic case 

The results in this section are new, but follow straightforwardly from previous work in Aoki 
et al and Aharony et al [40, 34]. Again we expand our fields, 



- oo 

A ab (r) = A«5 ab + Y, 

V Z7T 

v m=— oo 

v m=— oo 

- oo 



m=— oo 

but note that now with periodic boundary conditions we must include fermion zero modes 
As before we take 8A^, 8x\ °^ b , 8^^ ab to vanish when a = b, and the SU(N) colour 
symmetry implies that the sums ^2 a A a ,^2 a Xi and also X^a£a vanish. Using the notation 
A£ ab = £ a — £ 6 for the fermion zero modes we may now write the action to quadratic order 
in the perturbations, obtaining, 

squad _ ^ V V (S^ M> (m) ^ ( M ^ ab \ ( ^ ) (B 22) 

0<6m=— oo \ ^>Q u P / \ u *a,ab/ 

We see now that there are off-diagonal terms that mix the fermion and boson fluctuations, 
which are coupled together by the presence of the fermion zero modes. These can be removed 
by the following transformation, 

= ^i? - ^7-^^°^)^^)^ (B.23) 

which we note is a simple shift in the fermion fluctuation, and hence does not alter the path 
integral measure. Now we find the quadratic fluctuation action is diagonal, 

XT 00 f T\A m ) ab i n(m)ab „ \ / xff.(m) N 



a<bm=-oo \ P J \ ab 



(B.24) 
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where, 

S t» )ah = ; (1, (A^ ab a^) (, a jM at )(/Aa . (B.25) 

It follows from 5&^ b = and D^ ab = —D^ m " >ba that the action actually projects 

only onto the antisymmetric component S^J^ ab = \ (^sffi ab — sffi ab ^J . This antisymmetric 
part can be written as, 

T-,(m)a6 

s (m)ab = „ va p Lfa ( A ga6^ A ^a6) ( B . 26 ) 



We may now write the quadratic action as, 



m=-oo \a<6 a, 6 

P = 5% (D^D, - S^D 2 a + e^apDaJp) 5<f>„ (B.27) 

where, 

«/ a = -2i^p(Ae> a A0 (B.28) 

and in these last two equations we have suppressed the Fourier index m, and colour indices 
a, 6 for clarity. Again the bosonic operator in this action has a zero eigenvalue corresponding 
to the gauge freedom. Correspondingly we see, 

( M Mo6 + S (rn)^ D (rn)ab = q ( R2 g) 

Performing the same gauge fixing as before we have, 

S quad = ^_ j Y p9j ' + Yl iSA ° PD P SA ] - £ hi | det D T | 2 

m=— oo \a<b a,b I a<b 

P9 f- = (DiDj - SijDl + e lja pD a J p ) 5<S>j (B.30) 

and now integrating over the bosons yields the determinant of the 3x3 matrix in ij above, 
giving the elegant result 

OO j.j oo 

S quad = £ £ln(( J D 2 ) 2 + ( J D 2 )(J 2 )-( J D M J^) 2 ) + — Y, iSA ° PD ? SA ( R31 ) 

m=— oo a<b m=— oo a,fe 

where we again see that the SO (4) spacetime symmetry for the m = constant modes is 
restored after the gauge fixing had broken it. Performing the fermion integration then yields 
the 1-loop action, 

S quad = E £ln(l + ^((I^(J 2 )-(ZVn 2 )) (B.32) 
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and we immediately see that if it where not for the fermion zero modes contributing to the 
presence of J^, the action would vanish. Hence there has been almost total cancellation 
between the boson and fermion determinants. Now since ~ A£cj^A£ and each spinor £ Q 
has only 2 components, the action can contain at most quadratic terms in J^. Hence we may 
expand the logarithm above to give, 



gquad 



+ 



E E 

n=—oo a<b 
oo 

E E 

i=— oo a<b 



(jD (m)a6 )2 



n (m)ab n (m)a& 
4 \ U a u fi 



ab \ 



24 



ab a tob , 



(D { ; n)ab 



t-ab\ 



(B.33) 



where we have written the expression out fully. The partition function is then given by the 
functional integral over the bosonic and fermionic zero modes of this action, 

Z ~ / dAdX l d^e- S9uad 



I 



dAdX.dtdi 



n n 

i=—oo a<b 
( 



1 + 



24 



D 



{m)ab^ 2 



a<b V m=-c 



24 



3(AffAffA£fA£ 2 ab ) 



r(AefAg 6 AC? 6 AC 2 a6 ) 



(jD (m)ab )2 



(B.34) 



where we will not keep track of the overall normalization of Z. It is implicit in the measure 
that the gauge group is SU(N) and hence there is the constraint ^2 a £,aiY^ a £>a a = 0, an< ^ 
this is important in giving the form of the expression below. For example, if we had instead 
taken the gauge group U(N), the partition function would vanish as the trace of the adjoint 
fermion matrices that are constant on the circle decouple and since give fermion zero modes 
at 1-loop. After integrating these, the partition function will vanish. 

Following Aoki et al [40], we may perform the integral over the fermionic zero modes 
to give an effective 1-loop action for the bosonic moduli. In their analysis they express the 
answer in terms of sums over maximal trees with various valences as they consider 4,8 and 
16 supercharge matrix models. However as they note, the case of 4 supercharges is somewhat 
simpler than the others, and the 'tree' machinery is somewhat redundant. The answer can 
in fact simply be written as, 



dAdXi ^ M aia2 M a2CL3 . . . M aN - 2aN - 1 M aN - iaN 

(ai,a2,...,ajv)SP 

where P is the set of permutations of (1, 2, . . . , N) and, 
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This gives rise to the effective action given in the main text. We note that the infinite sum 
over Fourier modes is finite, and can simply be evaluated explicitly, although the expression 
obtained is unilluminating. 



C. Simulation details 

Both naive and supersymmetric lattice actions take the form 

S = k(S b (X) + S f @,^,X)) (C.l) 

where 

M 3 1 3 

x=0 i i>j 

has the same form in both cases except for a different definition of the covariant derivative. 
Notice that the scalar fields in the above action are expanded on the traceless antihermitian 
matrix basis of SU(N). The lattice coordinates x are equally distributed with spacing a on a 
circle of length R = Ma and R can be identified with ^ in the case of non-zero temperature. 
The fermion action takes the generic form 

S F = Y,^)M{X) x ^{y) (C.3) 

where the fermion operator A4(X) depends on the discretization. An immediate question 
arises as to how to scale k as the lattice spacing is reduced and the continuum limit ap- 
proached. Clearly the relevant dimensionless parameter in the continuum is Ai? 3 where the 
't Hooft coupling A = gp hys N is used to access the large N limit. Equating the inverse of this 
parameter to the corresponding lattice quantity yields the needed scaling of k 

NM 3 NM 3 . 

k = Uv=Hp- (c - 4) 

To simulate this system we must first integrate out the Grassmann fields yielding either 
det(M) or det5(M) for naive or supersymmetric actions. In practice this is accomplished by 
introducing complex commuting pseudofermion fields F with the same quantum numbers as 
the fermions ip and modifying the action as 

S = K^S B + F\M j! M)- p F S j (C.5) 

where p = \ or p = \ in the two cases. Notice that this requires the determinants be real 
and positive definite. 

Thus we now require an algorithm which can efficiently handle fields coupled through a 
fractional inverse power of the fermion operator - realized as a P x P matrix with dimension 
P = 2(N 2 — 1)M in the naive discretization (and twice this in the supersymmetric case). 
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To proceed further requires that we use an approximation for the fractional power. Most 
effective is a partial fraction realization of a rational function approximation in some interval 

Q 

x p ~a + Y e < x < 1 (C.6) 

f-f x + bi 

The optimal coefficients {ai,b a } can be determined offline using the remez algorithm which 
seeks to minimise the absolute value of the relative error for fixed (Q, e). In practice we have 
used approximations with Q = 10 — 15 and intervals ranging from e = 10~ 12 — 10 which 
conservatively covers the range needed and yields relative errors O(10 -4 — 10~ 8 ). The latter 
systematic error is far below the statistical errors of the Monte Carlo calculation and can thus 
be ignored. The action we have simulated thus takes the form 

S = k (sbM + apFtF + g aart-j—^Fj (C.7) 

which resembles the contribution of a number of doublets of degenerate fermions each with 
different mass parameters. 

This is still a non-local action and to simulate it efficiently requires the use of an auxiliary 
classical dynamics. To be precise one replaces the original partition function Z = J Dqe~ s ^ 
depending on a generic set of fields denoted by {q} by another comprising {q,p} with partition 
function Z' = f DqDpe~ H where 

H = s + E l p2 (c - 8) 

Clearly expectation values derived from Z' and Z are identical. Furthermore (approximate) 
Hamiltonian evolution in this phase space can be used to generate a series of global moves 
on {q,p} which, when subjected to the usual metropolis test, will generate the canonical 
ensemble needed to simulate Z' ? The resulting algorithm is termed Hybrid Monte Carlo 
[49]. 

In detail one starts from some initial set of coordinates q, draws new momenta p from a 
gaussian distribution and then evolves the fields {p, q} according to Hamilton's equations 

t ^ * <"> 

dp _ _as _ 
dt dq 

In practice a leapfrog integration with finite timestep dt is used to advance the fields along a 
classical trajectory of length r 

Pra+l/2 = Pn + (CIO) 
q n+ i = q n + Atp n+1/2 

Pn+1 — Pn+l/2 + "TT/ra+l 

3 This requires that the finite time step classical evolution be reversible and area preserving which is true 
for the leapfrog integrator used here 
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In practice we have used this update ut times corresponding to trajectory lengths r = 
rtyAt = 0.05 — 1. After this we expect the lattice Hamiltonian to be conserved to 0(Ai 2 ). 
We can then remove this At dependence by treating the final configuration as a potential 
global update in a metropolis simulation and accept the new configuration with probability 
e~^ H . Subsequently the momenta are refreshed and a new trajectory commenced. 

In our case q = {X, Ui,ip} and the major inputs to this algorithm are the forces 
J^- and Of these the most computationally expensive are the pseudofermion terms. For 

r) 

example the force term —ajf- is given by 



The final trick required to render this approach computationally feasible is to utilize a multi- 
mass solver to solve this set of Q equations iteratively and with a computational cost similar 
to the case when Q = 1. We have implemented a multimass CG-solver for our work [50]. 
Implementing the rational approximation for the fractional fermion operator in conjunction 
with the multimass solver in the way described is termed the Rational Hybrid Monte Carlo 
algorithm [51]. 

D. Continuum limits 

In this section we show the bosonic action observable Sb and the Polyakov loop observable 
P for increasing lattice sizes using the naive implementation for the quenched (figure 10), 
periodic (figure 11) and thermal theories (figure 12). We see that in all cases, over the range 
of (3 of interest, 5 lattice points already appears to give results that are close to the continuum 
limit. Note that for the periodic theory, the continuum limit of the bosonic action is discussed 
in the main text in the results section. 

We have observed that the thermal theory exhibits large lattice artifacts for small values 
of N and the number of lattice points M. This is true both for the naive and supersymmetric 
actions. To illustrate this figure 13 shows a plot of the Monte Carlo history of the maximal 
fermion eigenvalue A max for the SU(2) theory at P 3 = 10.0 obtained using the supersymmetric 
implementation. The plot includes data for lattice sizes M = 5, 10,20. The data is shown as 
a function of physical time measured in units of r = 1 RHMC trajectories. Data is plotted 
every 100 units of time with short time fluctuations being removed by plotting a running 
average obtained using a temporal window of 10 units of time. This fermion eigenvalue is 
strongly correlated with the rms value of the scalar eigenvalues. Very large fluctuations are 
seen with extremely large autocorrelation times for M = 5 lattice points. These fluctuations 
are not manifest in the bosonic action Sb and appear to derive from large fluctuations of the 




(C.ll) 



i=l 



where x % 1S the solution of the auxiliary problem 



(M*M + h) X l = F 



(C.12) 
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Figure 10: Plots of quenched theory action Sb, Polyakov loop P and scalar width W for N = 5 with 
varying lattice points M = 5,8, 12. 



scalars in the classical moduli space. These motions are strongly suppressed in the periodic 
theory on account of the presence of superpartner fermion zero modes. Of course as the 
number of lattice sites is increased the masses of these would be zero modes are lowered and 
once again they act so as to inhibit fluctuations of the scalars in the zero mode directions. We 
see this in the data which shows a marked reduction in the amplitudes of these fluctuations 
as the number of lattice points M increases. Of course the dimension of the moduli space 
varies like ./V and so we would also expect these effects to be suppressed at large N which is 
also observed. 

The problems seem most acute with the supersymmetric action and so for data presented 
in the main text we have concentrated on using the naive implementation for the bulk of our 
thermal runs. In this case to avoid these strong artifacts we require N > 5, M > 5 for the 
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Figure 11: Plots of periodic theory Polyakov loop P and scalar width W for N = 3 with varying 
lattice points M = 5, 8, 12. 



range of f3 we are studying, /3 3 < 10.0, and as we see in the figure 12 we obtain a good 
continuum limit. 

We conclude with figure 14 which shows data from the naive implementation of the 
fermion action Sf for the naive theory for representative values of N and 5 lattice points. 
Scaling arguments imply this quantity should equal unity and we see the data confirms this 
to within statistical errors, providing a non-trivial check of the RHMC routines. A similar 
check was performed on the supersymmetric implementation. 
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